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VIBRATION CHARACTERISTICS OF A LARGE WIND TURBINE 


TOWER ON NON-RIGID FOUNDATIONS 

by Suey T. Yee, Lewis Research Center, Cleveland, Ohio 
and Tse-Yung P. Chang, R. J. Scavuzzo, 

David H. Tiininerman , and John W. Fenton » 

Karl R. Rohrer Associates, Inc., Akron, Ohio 

SUMMARY 

A study was performed to determine the vibration characteris- 
tics of the Mod“OA wind turbine supported by a non-rigid foundation. 
The foundation considered is a square f coting with 34 feet in hori- 
zontal dimensions and 4 feet in thickness. Only the rocking action 
of the foundation was considered in the study. 

To establish a reference point, the free vibration of the wind 
turbine sitting on a fixed base was obtained from the dynamic analy- 
sis using the NASTRAN Computer program. Then, a simple model, con- 
sisting of three masses and springs representing the tower, machinery 
and blades, and a torsional spring representing the foundation-soil, 
u was used. The natural frequencies of the wind turbine were obtained 

for a foundation on soils with various rigidities. From the dynamic 
analyses, it was found that the Influence of foundation rotation on 
the fundamental frequency is quite significant for soils with elastic 
moduli less than 5000 psi (e.g., cohesive soil or loose sand) and the 
reduction in natural frequency can be greater than 20 percent. How- 
ever, for stiff er soils, such as well graded, dense granular materials 
or bedrock, the effect of the foundation-soil is not significant and 
therefore a rigid base can be assumed for dynamic analysis. 

In addition, the foundation effect with different tower heights 
was also studied. The analysis shows that the effect of soil flexi- 
bility on the vibration characteristics of the wind turbine is more 
pronounced if the height of the tower is increased. 
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lOTRODUCTlON 

With the declining supplies of oil end natural gas, alternative 
energy sources and energy conversion systems are being sought and 
an extensive energy research program is being developed under the direction of 
Energy Research and Development Administration (ERDA) . One of the 
many energy sources being examined is wind energy which has been 
utilized by mankind for centuries. The objective of this program 
is to assess the technology requiren«nts for constructing large wind 
turbine systems and ultimately to develop a wind energy system which 
is costocompetitive. In connection with ERDA's effort, NASA Lewis 
Research Center has designed and constructed a lOO-kilowatt wind 
turbine (called Hod*0) at its Plum Brook Station near Sandusky, 

Ohio. Larger wind turbines are being designed and will be built in 
other parts of the country where high wind zones are located. 

Xn designing the Mod~0 Wind Turbine, the tower was assumed to 
be resting on a rigid foundation; therefore, no interaction effect 
between the tower and foundation was considered. This assumption 
was justified by the fact that the Mod-0 tower is supported by 
fairly rigid concrete caissons sitting on hard shale. For wind 
turbines located on sites yet to be determined, the assumption of 
a rigid foundation may not apply. As a result, a study was initiated 
to determine the vibration characteristics of the wind turbine af- 
fected by non-rigid foundations. 

The general configuration of the proposed larger wind turbine 
(called )tod-OA) , similar to the one built at Plum Brook, consists of 
a tower; a nacelle housing the alternator, gear box and other 
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machinery; and two large rotor blades. The wind turbine is designed 
to produce 200 kilowatts of electric power in a 22*mph (hub-height) 
wind at a rotor speed of 40 r£»B. In the design of such a large 
structure f it is necessary to consider the dynamic forces caused 
by strong wind and forces induced by the rotation of the blades. 

Xn addition, the tower must be designed with sufficient rigidity 
such that its fundamental natural frequency is well above the 
critical driving frequency of the blades (1.33 Hs for a two blade 
rotor at 40 rpm) . Free vibration analyses of the blades fl) and 
the wind turbine as a whole 12) have been conducted for the Plum 
Brook system. A test program (31 was also performed to verify the 
analyses and it was concluded that the results obtained compare 
very closely with the test data. 

Included in this report is a study on the vibration characteris- 
tics of the Mod-OA wind turbine with a square footing on various 
soil conditions. A simple model consisting of three- lumped masses 
and springs was used for the study and only the rocking action of 
the foundation was considered. In order to establish a reference 
point, the free vibration of the wind turbine on a fixed base was 
obtained. Then, the description of the simple model and calculation 
of foundation spring constants were outlined. Finally, the natural 
frequencies of the wind turbine were obtained for various soil 
rigidities and the significance of the numerical results was 
discussed. 
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DESCRIPTION OF TOWER 

Tlui lto<l*OA wind turbine ccmsiete of a tower# a nacelle and 
t%fO rou>r blades# shown in Fig. 1# similar to the one Imiilt at 
Plim Brook* The tower is a 93-ft. tall welded trues 
constrvKited from steel pipe# angle and tee sections. The 
sectional view of the tower is a square shape and it varies from 
30 f^. wide at the base to 6.7 ft. wide at the top. Elevation 
of the tower and the sizes of sc»ne major m^nbers are shown in 
Fig. 2. 


The tower is anchored to a reinforced concrete slab foundation 
which is 34 feet square in horizontal directions and 4 feet in 
thickness. The dimensions of the foundation were based solely on 
the wind and dead loads, h nacelle# idiich contains the alternator, 
gear box# and other machinery# is attach^ to the top of the tower. 
Connected to the nacelle is a two blade rotor which is 125 feet in 
diazttter. The total weight of the nacelle and blades is 
approximately 44 #400 lbs. 


5 


FREE VIBRATION OF WIND TURBINE ON A FIXED BASE 

To oftabliBh a rofarence case for later diaeuasion, free 
vibration analyais of the wind turbine resting on a fixed base 
was performed 1;^ using the HASTRAN computer program [4]. The 
c^puter model consists of all the tower m&abers (represented 
by beam elements) and the masses due to the nacelle and blades 
which are evenly distributed at four corner nodes at the top of 
the tower. A total of 235 elements and 98 nodes were used. From 
the counter analysis, the first five fundamental frequ^mcies are 
tabulated in Table 1. 

The mode shapes for modes number 1, 2, 4, and 5 are plotted 
in Fig. 3; mode number 3, representing the vibration of local 
members, is not shown. The first and second modes represent the 
cantilever beam action of the tower in two independent x- and 
directions, respectively. The small difference in frequency of 
these two modes is due to the slight unsynenetric arrangement of 
the bracing number located at the second platform of the tower. 

The fourth mode shape represents the torsional action of the tower 
and the fifth is a higher bending mode. 

This analysis shows that the fundas^ntal vibrational mode is 
the cantilever beam action, particularly due to the heavy mass 
located at the top of the tower. Analyses were also conducted to 
investigate the ii^>ortance of rotary inertia of the blades on the 
fundan^ntal frequencies of the system. Numerical results indicated 
that this effect is very small, and therefore it can be ignored in 
the natural frequency analysis. 
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A SIHPLg HODEL 

Since the main objective of thie study is to detemine the 
effect of a non-rigid foundation on the vibration characteristics 
of the wind turbine, it is sufficient to use a sio^Xified model 
consisting of lumped masses and springs for the intended study. 
Only the rocking action of the foundation is considered. Based on 
the free vibration analysis of the tov%r on a fixed base, the 
fundamental vibrational mode is dus to the bending action. 
Therefore, it was decided that a model, consisting of 3-masses and 
springs representing the superstructure (the tower, nacelle and 
blades) and an effective spring representing the fouiulation-soil 
(shown in Fig. 4) be used, ^cording to this idealization, the 
flexibility and mass matrices are derived as follows. 

Consider in Fig. 4 that a horizontal force F^ is applied on 
tbm mass m^ causing a deflection Uj at fflj, then 

u^ ■ uj + uj » j - 1»2, or 3 (1) 

^ere u* represents the relative deflection at j before the rigid 
body motion; and Uj, the rigid body displacement due to the 
rotation of the tower. From the deflection analysis, one can 
easily find 


“j ■ 

(2) 

h h 

“J ■ 'i 

(3) 


X 


where fjj is the fie: i^llity coefficient of the tower structure on 
a fixed base; h^, elevation of the i^ mass; and X, the rotational 



•ti££n«ss o£ tha foundation. From Eq. (1)> on« can dofine the 
flexibility of the system as 
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C4) 


(5) 


The above equation consists of t\to terms: the first term is the 

fixed-base flexibility hnd the second term is due to the rotation 
of the foundation. The terms involving f® j were found from the 
deflection analyses of the actual tower using NASTRAH by applying 
unit loads at locations 1« 2, and 3< respectively. For the model 
TOnsidered, the mass points were selected at elevations 93-£t 
(location 1) , 68-ft (location 2 } , and 38-ft (location 3) . In 
this manner, the terms for f*j and f Jj are given by 




1 

ns 


972 456 

90 

285 

49 

symmetric 

14.5 



( 6 ) 


and 


Iff^l 


1 

K 


12.45 9.10 

5.10 

€.65 

3.72 

syBSMtric 

2.08 


in 


IF 


m 


In the above equation, determination of the foundation stiffness 
K will be discussed in the next section. 


The mass matrix of the model was determined by distributing 
the Blasses of the tower lumbers bet%^en two adjacent elevations 
equally to the respective controlling stations (i.e. 1* 2, or 3) 


Xa additioiif tha nassat of tha mchinary. and bladaa nara lua^ad 
at tha t<^# ioeatiofi 1 * ^arafora, the maas matrix waa found to 

ba 



N 


«#hara ■ 50.10 kipa 

^2 ■ 10.37 kipa 

W 3 - 16.27 kipa 

■ gravitational accalaration 


<•) 
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EFFECTIVE SPRING CONSTANT FOR FOUNDATION 

O^nsider a circular foundation of radius r^ resting on an 
elastic half space as shown in Fig. 5. The foundation is 
subjected to a constant moment T^i the corresponding amplitude 
of dynamic rotation may be evaluated from [5] 


A 
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V, My 


(9) 


where fg is the static rotation of the foundation and is the 
dynamic magnification factor. For a circular foundation of radius 
r^, the rotation is related to the moment by 


¥ 


S 


3(l-v) 

8 



( 10 ) 


where v is the Poisson's ratio and G, the shear modulus of soil. 
For foundations which are not circular « an effective radius based 
on equal areas can be found by 



Thus for a square foundation « Eq« (10) becomes 


u 


2.09 (1-v) 
G 



( 12 ) 


However, for a square foundation, a more exact expression is given 
by I5J 

1.92 (1-v) 


¥_ • 


G 


(13) 


Sq. (13) has bssn used to svaluata tha spring constant of tha 
foundation. 
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According to refaranca (5) tha dynamic magnification factor, 
Ny in Eq (9) is a function of tha mass ratio givan by 


B 


t 


3(l-v) 

8 



(14) 


and of tha dimansionlass frequency a^ defined by the equation 



(15) 


Where denotes the mass momnt of inertia of the foundation; p, 
tha mass density; tha forcing frequency; and the shear wave 
velocity of the soil media. Listing the equivalent radius based on 
Eq. (11), both By in Eq. (14) and ao in Eq. (15) can be easily 
determined and tha magnification factor My calculated. It was 
found that My is near unity * for tha frequency range of 
interest, i.e. f « 2.40 Hz. Therefore, the dynamic effect can 
be neglected, ftom Eq. (13) tha effective rotational spring 
constant K is 


K * 0.51 



( 16 ) 


where 


6 


E 


2(lw) 


E « Young's modulus 
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Easad on tho above eqiiation, the effective spring constant of the 
foundation can be readily calculated if the Young's modulus E and 
Poisson's ration v of soil is given. In the present investigation, 
a range of the Young's modulus for typical soils^ is considered as 
shown in Table 2. A constant Poisson's ratio of 0.3 is used in all 
calculations. Typical values of the effective foundation spring- 
constants vs. various soil nK>duli and shear wave velocities are 
tabulated in Table 3. 

j 


* Determination of My can be found in Appendix A. 

^ Determination of the foundation-soil parameters is outlined in 
l^)pendix B. 


MATURAL FREQUENCY ANALYSIS 


Referring to the simple model proposed# the equations of 
notion for the system can be written as 


If] [m] (u) + (u) ■ (0) 


( 17 ) 


idiere u denotes the absolute displacement vector# i.e. 

T 

(u) ■ (uj^, U 2 » U 3 > 


(18) 


and u is the corresponding acceleration vector. The matrices 
f and M were defined in Eqs. (4) and (8) respectively. 

The natural frequencies can be obtained from the solution to 
the following linear equation 


»^[f) CmJu) * (u) 


(19) 


where u is the angular frequency and u is the corresponding 
vibrational mode shape. The Iteration method [6] was employed 
to calculate the natural frequencies for different shear wave 
velocities and the results are tabulated in Table 4. As seen in 
the table# the first fundamental frequency obtained from the 
single model on a fixed base is 2.43 Hz which is very close to 
the value from the NASTRAN dynamic analysis of the tower, i.e. 
2.40 HZ. 

Define a frequency ratio for the first mode 
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where 

■ First natural frequency of the tower on a 
non-rigid base (spring) 

f^^ ■ First natural frequency of the tower on a 
fixed base. 

Thenf the frequency ratios for various soil rigidities are 
plotted against the shear wave velocities and elastic moduli in 
Figs. 6 and 7 respectively. As shown in these two figures, the 
change in natural frequency of the wind turbine on a non-rigid 
foundation varies widely with the stiffness of the soil. For 
dense, well graded granular materials or sound bedrock for which 
the elastic modulus exceeds 10,000 psi, the effect due to the 
foundation is very small, and therefore can be ignored in the 
vibration analysis of thw eind turbine. For medium to dense sand 
(E ■ 7000-12000 psi) , the percentage change in fundamental 
frequency (as con^ared with the fixed-base frequency) ranges 
between 4 and 18%. For cohesive soil or loose sand where the 
modulus is less than 5000 psi, the foundation effect is quite 
Important and the change in natural frequency can be greater than 


20 %. 
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APPROXIMATE ANALYSIS OF FOUNDATION EFFEC TS 

Tha effects of tower flexibility and foundation rotation on 

the fundamental frequency can be evaluated separately in an 

approximate manner. By assuming that the tower is fixed at the 

foundation, the fundamental frequency caused by tower flexibility 

can be calculated. This condition is obtained when v ■ • and 

s 

results are listed in Table 4. For the tower investigated fj^^ <■ 

2.40 Hz. 

The effect of tower foundation rotation only on the fundan^ntal 
frequency of the structure can be evaluated assuming that the tower 
is rigid. For this case the frequency of the system can be 
calculated as follows 



where 


K ■ effective foundation rotational stiffness 
(Table 3) 

I - Zmiti 

For this tower a soil shear velocity of V * 250 ft/sec, and 

s 

frequency as determined from Eq. (19) is f^^^ ■ 


3.13 Hz. 



These two effects, tower flexibility and tower rotation, 
can be combined using Dunkerly's formula (Ref. 6) to estimate 
the fundamental frequency of the system. 

JL - J- + 

(^l) ^ (^Ir^ (^Is) ^ 

For this case f = 1.92 which is close to the more exact value 
listed in Table 4. 

The advantage of this approximate method is that effects of 
foundation rotation on the tower fundamental frequency at a 
particular site can be easily estimated prior to a more complete 
analysis. These equations show that as the tower becomes taller 
or stiffer the effect of the soil is more pronounced, if other 
factors remain constant. 


CONCLUSION 


Vibration characteristics of a 200-kilowatt wind turbine 
system resting on rigid and non-rigid foundations were studied. 
Ba&;ed on the study, the following conclusions can be drawn: 

(1) For a rigid foundation, the fundamental mode of the 
tower is 2.4 Hz and is caused basically by cantilever 
bending action of the tower. 

(2) The influence of the foundation rotation on the 
fundamental frequency is quite significant for soils 
with elastic moduli less than 5000 psi. 

(3) For stiffer soil, such as dense, well-graded granular 
materials or sound bedrock, the effect of foundation 
rotation is not significant and the effect can be 
ignored. 

For other wind turbine tower designs in which the height of the 
tower is increased while other parameters are unchanged, the 
approximate analysis in the previous section indicates that the 
vibration characteristics of the system become more sensitive to 
changes in soil flexibility. Furthermore, if the tower design is 
altered to increase its stiffness, the influence of the soil 
flexibility may dominate the dynamic response of the system. 


17 


APPENDIX A - DETERMINATION OP DYNAMIC MAGNIFICATION FACTOR 


The dynamic magnification factor M^ in Eq. (9) is. a 
function of the mass ratio and a dimensionless frequency a^ 
defined by 

3(l-v) l^ 

- £-2- (A-1) 

8 P 


ur 

o 


a 


o 


Since the fundamental frequency of the wind turbine is 
« 2.43 Hz. 

and the width and the length of the foundation are 
L - d = 34’ 

Then u> * » 15.27 rad/sec 

r^- 230 in - 19.17 ft 

The mass ratio, B^, can be estimated as follows 



(A-2) 


(A-3) 


where 

W » weight of foundation 
t » length 

Since the foundation is to be constructed on concrete aE^roximately 
4 feet thick and 34 feet square 
W - Y.V_ 

^ 2 3 

W - (150) j^3 (34)^ (4) ft^ 

W - 693,600 lbs 


(A-4) 
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Thus 


If O 

I. ■ : < > * 2.08 X 10* Ib-ft-sac 

’ « IT" 


According to Eq. (A-1) 


(A-5) 




(3) (0,7) (2.08) (10 ) 


( 8 )( 


100 

32.2 


230 5 

) ( — ) 
12 


0.068 


(A-6) 


Using Figure 7-15 of Reference C 1 1 # it can be seen that 
is unity for the frequency range of interest. Therefore » 
this effect can be neglected. 
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APPENDIX B - DETERMINATION OF FOUNDATION SOIL PARAMETERS 


Typical analysis techniques for determining the response of 
soil-foundation systems subjected to dynamic loadings consider 
the soil as son« type of equivalent elastic system. Since soil 
is not a linearly elastic material, an approximation must be 

made of the elastic modulus, E, or the shear modulus, G, and 

% 

Poisson's ratio, v, which will produce a calculated response 
within reasonable accuracy. 

1. Poisson's Ratio 

Poisson's ratio is difficult to determine accurately 
for soils and can be assumed for most practical calculations 
based on knowledge of the boil type. Typical values are as 
follows : 

Saturated Clay v - 0.50 

Unsaturated Clay and Clay 

with Sand and Silt v ■ 0.30 to 0.40 

Granular Soils v ■ 0.30 to 0.35 

A value of 0.3 has been used in Section 3.0. 

2. Elastic Modulus 

The elastic modulus and shear modulus vary over a wide 
range depending on soil type, density and confining pressure. 
Typical values of the elastic modulus, E, are listed in Table 2. 
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3. Laboratory Testing and Field Testing 

Determination of actual values of E or 6 can be made 
from emperical relationships, frcMs laboratory testing or 
preferably from in situ measurements. 

Empirical relationships have been developed by Richer t 
and Hardin [5] for determining the shear modulus, 6, for 
granular soils. 

For round^grained sands (e < 0.80) i 
2630(2. 17-e)^(3)®*® 

G - 2 (B-1) 

1+e 

For angular-grained materials: 

1230(2. 97-e)^(Oo) 

G - 

14-e 

where, G ■ shear modulus - psi 

0 ^ ■ average effective confining pressure 
e • void ratio 

Laboratory testing is normally performed by the resonant 
column method. Undisturbed samples of cohesive soils can 
readily be tested by this method. However, undisturbed 
san^les of granular soils cannot easily be obtained and since G 
depends on void ratio and confining pressure, determining G 
from laboratory testing is usually not practical. 


(B-2) 
- psi 
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Khenever possible « In situ evaluation of the shear 
modulus should be carried out. This Involves determining 
the average shear wave velocity# for the foundation 
supporting soil. The shear modulus can then be determined 
from 

G - (B-3) 

where p equals the mass density of the soil. 

The shear wave velocity# V^# can be determined either 
by direct measuren«nt using seismic survey methods or preferably 
from steady-state vibration methods. The steady-state 
vibration method measures Rayleigh wave velocities which for 
practical engineering purposes can be assumed to equal shear 
wave velocities. By varying the frequency of the wave 
generating vibrator# the average shear modulus can be 
obtained for different depths Into the soil for use In 
analyzing possible soil- foundation systems. 
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TABLE I 

rUMDAMENTAL FREQUENCIES OF HIND TURBINE ON FIXED BASE 


Mod# 

Description 

Value 

X 

BeeiRwise (x direction) 

2.37 

2 

Beamvise (y direction) 

2.40 

3 

Local H«ftber 

6.67 

4 

Torsional 

7.29 

5 

Beamirise (y direction) 

11.07 


TABLE 2 

YOUNG'S MODULI OF TYPICAL SOILS 


Soil Type E(P«i) 


Very soft to soft clay 

50 - 

500 

Medium to hard clay 

500 - 

2,500 

Sandy clay 

4,000 - 

6,000 

Loose sand 

1,000 - 

4,000 

Dense sand 

7,000 - 

12,000 

Dense sand and gravel 

14,000 - 

28,000 



TABLE 3 
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EFFECTIVE ROTATIONAL SPRING CONSTANTS OF FOUNDATION 


V, 

G 

c 

K 

trt/s«c) 

(P*i) 

IZSil 

(Lb-In/Rad) 

250 

1,350 

3,500 

6.81 X 10^® 

$00 

5,390 

14,000 « 

27.19 X 10^® 

600 

7,760 

20,200 

39.15 X 10^® 

1000 

21,600 

56,100 

10.90 X 10^^ 

1500 

48,500 

126,200 

24.47 X 10^^ 

2000 

86,300 

224,300 

43.54 X 10^^ 


V. 

TABLE 4 

NATURAL FREQUENCIES 



Ft/t«c 


W 


(Kz) 

250 

1.89 

15.4 

% 

29.6 

500 

2.23 

18.5 


36.7 

1000 

2.35 

19.8 


46.9 . 

1500 

2.38 

20.1 


51.7 

2000 

2.39 

20.1 


54.4 


2.40 

20.3 


57.2 



AL 

iSoR QUW.V 


TOWER LEVELS (FT) 


Z 


93 

81 

68 

54 

38 
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FIG. 2 ELEVATION OF TOWER STRUCTURE 
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FIG. 3 MODE SHAPES (CONTINUED) 
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FIG. 3 MODE SHAPES (CONTINUED) 






FREQUENCY RATIO 


r-l 

I 


O 

H 

S 

>* 

u 

z 

M 

D 

a 

s 

(14 



DENSE 


LOOSE SAND DENSE SAND SAND & GRAVEL 



ELASTIC MODULUS, E (10^ psi) 


FIG. 7 FREQUENCY RATIO VS. ELASTIC MODULUS 




